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INTRODUCTION 


In  Technical  Report  Nos.  15,  "Ultrasonic  Wave  Reflection  at  Liquid- 
Solid  Interfaces,"  and  17,  "Studies  of  Linear  and  Nonlinear  Ultrasonic 
Phenomena,"  we  presented  summaries  of  our  contributions  to  certain 
subjects.  Technical  Report  No.  18  is  intended  to  expand  on  the  summary 
and  bring  it  up  to  date.  It  is  divided  into  two  parts. 

Part  I.  Schlieren  Studies  of  Ultrasonic  Waves 

The  Report  begins  with  the  description  of  a  unique  goniometer  designed 
by  members  of  the  Ultrasonic  Group  for  use  in  the  schlieren  system  for 
visualization  of  ultrasonic  waves  in  liquids.  By  using  the  properties 
of  parallelograms  we  were  able  to  produce  a  precision  goniometer  without 
use  of  precision  machine  shop  facilities.  The  second  paper  presents 
some  photographs  made  with  the  goniometer  in  the  schlieren  system  and 
shows  the  effect  of  a  layer  of  A^O-j  on  a  stainless  steel  reflector  of 
ultrasonic  waves  in  water.  The  leaky  Rayleigh  wave  excited  in  the 
A1203  layer  has  a  velocity  smaller  than  that  exicted  either  at  a  water- 
stainless  steel  interface  or  at  a  water-A^O^  interface. 

Part  II.  Nonlinear  Acoustics  of  Solids 

In  relatively  larqe  sinale  crystal  samples  (1  inch  in  diameter  and 
1  inch  long)  one  can  measure  such  things  as  "The  Nonlinearity  Parameters 
and  Third-Order  Elastic  Constants  of  Copper  between  300  and  3°  K"  as 
reported  in  Paper  No.  3.  The  fact  that  the  measurements  can  be  made 
to  low  temperatures  is  especially  important,  as  the  effect  of  thermal 
motion  of  the  atoms  is  ignored  in  many  theories.  This  means  that  they 
are  strictly  applicable  only  at  0°  K.  For  comparison  with  these  theories, 
then,  we  measure  to  the  lowest  readily  obtainable  temperature. 
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A  problem  encountered  in  the  study  of  the  nonlinear  properties  of 
solids  is  the  fact  that  oftentimes  it  is  difficult  to  grow  large  single 
crystals  of  interesting  substances.  Ordinarily  we  use  a  30  MHz  ultra¬ 
sonic  wave  of  finite  amplitude  to  determine  the  nonlinearity  parameters 
of  single  crystals  1  inch  in  diameter  and  1  inch  long.  The  amplitude 

of  the  second  harmonic,  which  must  be  measured  absolutely,  typically  is 

-2  0 

of  the  order  of  10  A  in  these  samples.  We  posed  for  ourselves  a 
question:  Given  our  desire  to  measure  nonlinearity  parameters,  what  is 
the  smallest  sample  one  can  measure  with  present  technique?  The  fourth 
paper,  "Measurement  of  Nonlinearity  Parameters  in  Small  Solid  Samples 
by  the  Harmonic  Generation  Technique,"  is  an  attempt  to  answer  the 
question. 

Another  question  of  fundamental  importance  to  nonlinear  acoustics 
of  solids  is  the  relationship  between  the  nonlinearity  parameter 
measured  acoustically  and  the  Griineisen  parameter  which  comes  from 
measurement  of  thermal  properties.  This  question  is  given  a  relatively 
simple,  and  almost  complete,  answer  in  the  fifth  paper  on  "Relationship 
between  Solid  Nonlinearity  Parameters  and  Thermodynamic  Griineisen 
Parameters."  This  paper  was  based  on  the  oral  presentation  given  at  the 
joint  meeting  of  the  Acoustical  Society  of  America  and  the  Acoustical 
Society  of  Japan.  This  was  an  especially  appropriate  audience  since  two 
of  the  authors  were  from  the  United  States  and  one  was  from  Japan. 

The  final  paper  in  this  Report,  "Quantum  Mechanical  Theory  of  Non¬ 
linear  Interaction  of  Ultrasonic  Waves,"  answers  in  part  another 
fundamental  question  of  nonlinear  acoustics.  Presumably  in  the  corres¬ 
pondence  limit  the  quantum  mechanical  description  of  phonon-phonon 
interaction  would  become  identical  to  the  description  (based  on  nonlinear 
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elasticity) of  the  scattering  of  one  acoustical  disturbance  by  another. 

But  this  assumption  was  hard  to  prove.  The  paper  provides  specific 
examples.  It  begins  with  the  general  quantum  mechanical  description 
of  phonon-phonon  interaction  and  specializes  the  description  to  that 
of  two  col  linear  phonons  of  frequency  v  which  interact  to  produce  a 
phonon  of  frequency  2v.  (This  comes  from  energy  conservation: 
hv  +  hv  =  2hv.)  By  maintaining  the  wave  description  (avoiding  quantiza¬ 
tion),  one  is  able  to  show  that  the  mathematical  result  is  identical 
to  that  previously  obtained  from  a  generalization  of  elasticity.  This 
is  true  also  in  the  description  of  third  harmonic  generation.  In  third 
harmonic  generation  one  is  able  to  show,  in  addition,  that  the  small 
term  in  the  third  harmonic  amplitude  which  contains  fourth-order  elastic 
constants  in  fact  comes  from  four-phonon  interactions  in  the  quantum 
mechanical  picture,  whereas  all  of  the  other  terms  (those  involving  second 
order  and  third-order  elastic  constants)  resulted  from  three-phonon 
interactions.  The  advantage  of  the  quantum  mechanical  approach  lies 
primarily  in  the  fact  that  the  path  from  the  general  description  to  the 
particular  application  is  explicitly  marked,  and  the  point  at  which  one 
makes  any  particular  simplifying  assumption  can  be  located  unambiguously. 
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A  UNIQUE  GONIOMETER  FOR  USE  IN  SCHLIEREN 
VISUALIZATION  OF  ULTRASONIC  WAVES 


M.  A.  Breazeale 
Department  of  Physics 
The  University  of  Tennessee 
Knoxville,  Tennessee  37916 


Introduction 

Alignment  problems  encountered  in  the  use  of  schlieren  systems  have 

been  recognized,  and  often  solved,  by  a  number  of  investigators.  For 

example.  Fig.  1  is  a  schlieren  photograph^  of  the  interaction  between  an 

? 

ultrasonic  beam  of  Gaussian  cross  section  and  a  leaky  wave  at  an  inter¬ 
face.  The  fact  that  the  reflected  beam  is  displaced  to  the  left,  rather 
than  to  the  right,  results  from  a  fine  periodic  structure  (grating)  at 
the  interface.  The  periodic  structure  shifts  the  phase  of  a  leaky  wave 
trapped  at  the  interface  by  180°,  and  causes  the  energy  flow  to  follow 
the  path  indicated  in  Fig.  2. 

The  photograph  of  Fig.  1  was  made  only  after  expenditure  of 
considerable  effort.  The  grating  periodicity  of  0.178  mm  reguired  not 
only  precision  alignment,  but  also  the  maintenance  of  the  precision 
alignment  as  the  incident  angle  was  changed.  Subseauently  a  movie  was 

3 

made'  to  show  in  detail  what  happens  as  the  incident  angle  is  changed. 
For  the  movie  the  point  of  contact  between  the  ultrasonic  beam  and  the 
reflecting  interface  had  to  remain  the  same  for  all  incident  angles. 

This  required  a  precision  goniometer  capable  of  rotating  the  transducer 
along  a  circular  path  centered  at  the  point  of  contact  between  the 
ultrasonic  beam  and  the  reflecting  interface. 
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Figure  1.  Schlieren  photograph  of  an  ultrasonic  beam  of  Gaussian  cross 
section  reflected  from  a  brass  grating  in  water. 


Figure  2.  Direction  of  energy  flow  in  ultrasonic  beam  shown  in  Figure 
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The  purpose  of  the  present  discussion  is  not  to  give  an  exhaustive 
description  of  the  phenomena  that  occur  when  an  ultrasonic  beam  is 
reflected  at  a  liquid-solid  interface.  Rather,  the  purpose  of  the 
present  discussion  is  to  describe  a  precision  goniometer  that  has  made 
an  experimental  study  of  such  phenomena  possible. 

Description  of  Goniometer 

The  design  of  the  goniometer  can  be  understood  by  observing 
geometrical  properties  of  parallelograms.  Suppose  a  parallelogram  were 
constructed  by  placing  bearings  at  the  four  corners.  The  parallelogram 
then  could  distort  as  shown  in  Fig.  3,  in  which  two  positions  are  drawn. 
The  property  of  interest  here  is  the  fact  that  any  point  on  the  top  edge 
of  the  parallelogram  traces  out  a  circle  centered  at  a  corresponding 
point  on  the  bottom  edge.  This  is  true  not  only  at  the  corner,  labelled 
A-j  and  also  for  anY  other  point  (e.g.,  B-j  and  B^),  as  indicated. 

The  centers  of  the  circles  are  and  Cg,  respectively. 

Let  us  now  remove  the  bottom  edge  of  the  parallelogram,  but  fix  the 
position  of  the  bearings  by  adding  an  auxiliary  support.  Further,  let 
us  form  another  parallelogram  by  adding  a  horizontal  piece,  as  indicated 
in  Fig.  4.  To  the  two  horizontal  portions  of  the  new  parallelogram,  let 
us  now  add  a  new  set  of  bearings  and  attach  a  new  vertical  member.  A 
transducer  attached  to  this  piece  has  the  capability  we  desire:  the 
possibility  to  rotate  about  a  point,  the  point  of  intersection  of  the 
two  dotted  lines.  A  reflecting  surface  placed  at  the  position  indicated 
can  be  studied  in  detail.  The  ultrasonic  beam  reflected  from  it  will 
reflect  from  exactly  the  same  point  on  the  interface  for  all  incident 
angles. 
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A  goniometer  using  the  principles  described  has  been  in  use  in  the 
Ultrasonics  Laboratory  at  The  University  of  Tennessee  for  several  years. 

A  photograph  of  the  goniometer  and  one  lens  of  the  schlieren  system  is 
shown  in  Fig.  5.  In  the  photograph  several  refinements  are  observable: 

1.  The  support  for  the  transducer  has  a  dog-leg  in  it.  This  avoids 
problems  that  would  arise  for  incident  angles  great  enough  that 
the  edge  of  the  water  tank  would  interfere  with  transducer 
motion. 

2.  A  scale  is  provided  on  one  of  the  supports  so  the  incident  anale 
can  be  measured  directly. 

3.  A  worm  gear  is  used  to  make  precision  adjustment  of  the 
incident  angle. 

4.  The  entire  system  is  counterbalanced  by  lead  bricks. 

Finally,  it  is  obvious  that  the  same  effect  could  be  accomplished 

by  use  of  a  circular  track  of  sufficiently  large  diameter.  Such  systems 
have  been  constructed  in  other  laboratories,  but  they  reguire  complicated 
and  expensive  machine  shop  work.  The  advantages  of  the  paralleloaram 
goniometer  are: 

1.  Simplicity  of  construction.  (One  only  needs  to  bore  the  holes 

for  the  bearings  at  the  correct  positions.  No  large  or  complicated 
machine  shop  facilities  are  needed.) 

2.  Precise  adjustment  possible  (no  stick-slip  in  the  movement). 

3.  Maintenance  of  precision  over  long  periods. 
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Figure  5 


Parallelogram  ultrasonic  goniometer  in  a  schlieren  system. 


Conclusion 


This  new  approach  to  an  old  problem  nas  produced  a  goniometer  that 
is  working  effectively  in  a  schlieren  system  having  an  aperture  of  8  inchec. 
The  principle  on  which  the  goniometer  is  based  can  be  used  with  even 
larger  systems  without  a  fundamental  increase  in  the  complication  of 
construction. 
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LEAKY  WAVE  GENERATION  AT  AL,Oj  LAYER  ON  '•'UNLESS  jTELl.  IN  WM I.R 


Laszlo  Adler  and  Daniel  L.  Butler 


Department  of  Physics 
The  University  of  Tennessee 
Knoxville,  Tennessee  37916 


Abstract 

When  a  finite  ultrasonic  beam  with  a  Gaussian 
amplitude  distribution  is  reflected  from  an  AI2O3 
layer  on  stainless  steel  in  water  the  reflected 
amplitude  field  distribution  indicates  leaky  wave 
generation.  The  angle  at  which  the  leaky  wave  is 
generated  defines  the  leaky  wave  velocity.  The 
leaky  wave  velocity  of  this  structure  is  measured 
as  a  function  of  kh  (where  k  is  the  wave  number 
and  h  is  the  layer  thickness).  The  effect  of  the 
presence  of  the  layer  on  the  reflected  amplitude 
distribution  is  discussed. 


1.  Introduction  and  Background 

The  reflection  of  ultrasonic  waves  at  a 
liquid-solid  interface  as  a  function  of  incident 
angle  is  a  basic  boundary  value  problem.  For 
infinite  plane  waves  one  solves  the  wave  equation 
with  appropriate  B.C.  Recently  interest  has  been 
focused  on  the  physically  realistic  problem  of 
the  reflection  of  a  finite  beam  of  some  well' 
defined  shape  from  the  liquid-solid  interface. 

Both  experimental  1  •  *  and  theoretical*^  analysis 
have  established  the  existence  of  the  so-called 
leaky  Rayleigh  waves  at  the  liquid-solid  boundary, 
The  existence  of  these  waves  is  easily  demonstrated 
by  a  Schlicrcn  picture  of  an  incident  Gaussian 
ultrasonic  beam  reflected  below,  at  and  above  the 
Rayleigh  angle.  Figure  l  shows  the  case  for  water- 
stainless  steel  interface.  The  middle  picture  is 
taken  at  the  Rayleigh  angle  about  30.5*.  The  angle 
of  incidence  is  such  that  the  refracted  wave  is 
coupled  along  the  interface  and  leaks  back  to  the 
liquid  as  it  propagates.  This  leaky  wave  together 
with  the  specularly  reflected  field  produces  the 
total  reflected  field. 

The  theoretical  analysis  of  this  finite  beam 
problem  wa»  carried  out  by  Bertoni  and  Tamir. '  In 
their  development  of  an  analytical  approximation  to 
the  reflection  integral,  Bertoni  and  Tamir  show 
that  it  divides  into  two  parts  and  can  be  written 
as 

wrefl(,’:)  ’  vo*x,J*  *  vj(*.zl  •  (I) 

vQ  (  * .  1 )  represents  a  specular  reflection,  which 
resembles  the  incident  beam  in  its  amplitude  dis¬ 
tribution  hut  is  shifted  180"  in  phase,  v  («,z)  is 
the  surface  component,  which  is  in  phase  with  the 
incident  heam  over  part  of  the  interface  and  out  of 
phase  over  the  remainder;  its  amplitude 


distribution  is,  in  general,  significantly 
different  from  that  of  the  incident  beam  and 
the  specular  component. 


Figure  1  -  Reflection  of  a  Gaussian  ultrasonic  Beam 
from  a  Water -St  a  in  1  ess  Steel  s  rtace  for 
Incident  Angles  (at  JS",  (hi  '  .S"  (Rayleigh 
Angle),  (c I  40’ 


1979  Ultrasonics  Symposium  Proceedings,  IEEE  Cat.  »79CII1482-9SU 


Laszlo  Adler  and  Daniel  L.  Butler 
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Tor  a  Gaussian  incident  bean,  us  used  in  the 
experiments  described  below,  Bertoni  and  Tam r-5 
obtained  an  analytical  approximation  valid  at  the 
interface.  In  order  to  compare  theory  with  experi¬ 
ment,  Brcazeale,  Adler,  and  Scott-  corrected  the 
approximation  for  points  in  the  liquid  halfspace. 
The  expressions  for  the  leaky  wave  field  compo¬ 
nents  are: 


vw  * 


j  (-  (xf/wr)  nk[xrsin9p*(zr-z0)cos9p]  i 

'  /iTw  cose 

r  P 

yVw  , 

l'lur,Ir)  *  -2V!  t  *r .  =r)  1 1  --ft-5-  exp(y2)erfc(y))  (3) 


where 


and 
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The  beam  halfwidth  w  is  measured  at  z,,  0_  is  the 
liquid-solid  equivalent  of  the  Rayleigh  ahgle,  and 
A  is  the  so-called  "Schoch  displacement."  A  was 
derived  by  Schoch  in  his  original  lateral  displace¬ 
ment  theory  and  was  shown  to  be  mathematically 
equivalent  to  a  surface  wave  decay  constant 
occurring  in  the  Bertoni -Tamir  approximation.  It 
is  a  complex  function  of  the  acoustic  velocities 
and  densities  of  the  interface  media.  The  finite 
beam  reflection  from  the  solid  layer-solid  in  water 
has  no  theoretical  treatment  at  present.  The 
experimental  investigation  of  this  problem  will  be 
presented  in  the  next  section. 


2.  Leaky  Wave  Generation  of  Solid  Layer- 
Solid  Interface  in  Water 


The  problem  presented  here  deals  with  the 
generation  of  leaky  waves  at  a  liquid-solid  layer- 
solid  interface.  The  problem  is  shown  on  Fig.  2. 
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Figure  2  -  Water-Solid  Layer-Solid  Interface. 

A  solid  layer  which  is  of  AI2O3  ceramic  layer  is 
sprayed  on  a  stainless  steel  block.  The  reason 


AlsOj  is  chosen  is  because  of  our  previous 
studies  of  leaky  wave  structure  on  both  water- 
Al20^  and  on  water  stainless  steel.-’  Ihe  main 
part  ot  the  experimental  arrangement  to  obtain 
quant .  at  1 ve  data  of  the  reflected  beam  profile  is 
a  specially  assigned  goniometer  shown  on  Iig.  3. 

A  2  MHz  Gaussian  transducer  sends  out  some  long 
pulses  of  20-30  usee.  The  receiver  is  scanned 
through  about  7  cm,  which  is  the  extent  of  the 
reflected  field. 


Figure  3  -  Ultrasonic  Goniometer. 


The  demonstration  of  the  existence  of  leakv 
waves  was  done  in  the  following  ways:  The  receiver 
and  transducer  arms  were  scanned  by  small  incre¬ 
ments  of  angle'  until  the  RF  waveforms  indicated 
the  phase  cancellation.  On  Fig.  1  the  RF  wave¬ 
form  reflected  from  the  A^Oi-stainless  steel  in 
water  is  shown.  The  angle  at  which  the  phenomena 
took  place  is  31.5“.  The  Rayleigh  leaky  velocity 
is  calculated  to  be  2.0.3  »  10s  cm/sec.  The  thick¬ 
ness  of  the  AljOj  layer  is  80c.  The  leaky  wave 
velocity  from  the  AljDj  laver-stainless  steel  was 
measured  as  a  function  of  kh.  On  Fig.  5  the  result 


Figure  4  -  Received  Signal  in  the  Ravlcigh  Angle 
Reflection  from  Water-AlaOj  Laver-Stainless 
Steel . 

is  plotted.  At  kh  »  0,  the  Rayleigh  leaky  velocity 
of  the  stainless  steel  is  obtained.  As  kh  -  -,  the 
Rayleigh  velocity  pf  the  sprayed  Al;0j  obtained  (in 
bulk!  is  2.18  x  10J  cm/soc.  which  is  significantly 
lower  than  the  velocity  of  the  AljOj  in  compact 
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form  (Vn  *  5.53  x  10^  cm/scc).  From  the  asymptotic 
value  of  the  leaky  wave  velocity  one  may  obtain  the 
shear  velocity  and  shear  modulus  as  well  of  the 
porous  material  AljO;.  Additional  knowledge  of  the 
wave  propagation  in  the  thin  layered  material  may 
he  obtained  from  the  measurements  of  the  ampli¬ 
tude  distribution  of  the  reflected  wave  at  the 


Figure  5  -  Variation  of  Leaky  Rayleigh  Velocity 
(VR)  for  Water-A^Oj  Layer-Stainless  Steel 
Interface. 


Rayleigh  angle.  On  Fig.  6  the  theoretical  values 
of  the  reflected  field  is  plotted  (solid  line). 

This  curve  was  calculated  from  the  modified 
Bertoni  theory  for  liquid-steel  interface  and  for 
the  parameters  used  in  the  experiment  (beam  width 
w  »  8  mm,  f  »  2  MHz,  transmitter- interface- 
receiver  total  distance  zR  »  400  ms).  The 
corresponding  experimental  points  measured  are  in 
good  agreement  with  the  theory.  The  effect  of 
the  80u  AI2O3  layer  on  the  experimental  data  is 
shown  by  the  points.  There  are  several  features 
of  these  latter  data  points  to  consider.  First, 
the  amplitude  of  the  first  peak  is  diminished 
compared  to  the  one  without  the  layer.  Second, 
the  null  point  which  indicates  the  phase  can¬ 
cellation  has  shifted  laterally,  indicating  that 
the  parameters  changed  and  at  phase  cancellation 
between  Vo  and  Vj  will  take  place  at  another 
point.  Since  at  that  point  Vo  *  Vj  equations  may 
be  inverted  to  obtain  additional  parameters  of  the 
interface  when  a  thin  layer  is  present  on  a  sub¬ 
strate.  Further  theoretical  work  is  required  to 
interpret  these  experimental  findings. 

3.  Conclusions 

The  problem  of  ultrasonic  leaky  wave 
propagation  in  thin  (80u)  ceramic  layer  (AI7O3)  on 
stainless  steel  immersed  in  water  has  been  Investi¬ 
gated.  It  appears  that  the  leaky  wave  velocity 
becomes  dispersive  with  the  presence  of  the  layer. 
The  layer  also  affects  the  amplitude  distribution 
of  the  reflected  field.  No  theoretical  work  is 
available  at  present  to  analyze  this  problem. 
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Figure  6  -  Amplitude  Distribution  of  a  Reflected 

Gaussian  Ultrasonic  Beam  at  the  Rayleigh  Angle 
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ABSTRACT 

The  ultrasonic  harmonic  generation  technique  has  been  used  to  extend 
measurement  of  the  nonlinearity  parameters  of  copper  to  3  °K.  Comparison 
of  these  data  and  combinations  of  truly  adiabatic  TOE  constants  with 
predictions  of  simplified  models  show  that  a  central  force,  nearest 
neighbor  model  accounts  reasonably  well  for  the  behavior  of  copper  in 
the  regions  of  45  °K  and  200  °K  and  less  well  at  other  regions.  The 
central  force,  nearest  neighbor  model  also  gives  a  good  qualitative 
explanation  for  the  temperature  dependence  of  the  combinations  of  TOE 
constants  that  are  measured  in  this  investigation. 


To  be  published  in  the  December  1980  issue  of  Journal  of  Applied  Physics. 
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I .  INTRODUCTION 

In  the  present  investigation  previously  established  techniques  are  used 
to  measure  combinations  of  third-order  elastic  constants  of  copper  from 
300  °K  to  3  °K.  The  technique  involves  the  measurement  of  the  distortion 
of  an  ultrasonic  wave  as  it  propagates  through  various  copper  single 
crystals. 

The  development  of  a  capacitive  detector  and  its  calibration 

permit  the  absolute  determination  of  the  amplitudes  of  finite  amplitude 

ultrasonic  waves. *  Later  refinements  made  possible  the  extension  of 

2 

these  measurements  to  lower  temperatures.  These  methods  have  been 
used  to  calculate  various  combinations  of  TOE  constants  at  low 
temperatures.  Peters,  Breazeale,  and  Pare3used  this  technique  to 
measure  combinations  of  TOE  constants  of  copper  to  77  °K.  Yost  and 
Breazeale'1  measured  combinations  of  TOE  constants  of  germanium  to  77  °K. 

Bains  and  Brcazcalc^extended  the  measurements  of  germanium  to  3  °K. 

Cantrell  and  Breazeale^measured  Cjjj  for  various  samples  of  fused 
silica  between  300  and  3  °K. 

Various  investigations  of  the  TOE  constants  of  copper  have  been 
7 

made.  Daniels  and  Smith  isolated  various  combinations  of  TOE  constants 

for  copper  by  measuring  the  pressure  derivatives  of  second-order 

8 

clastic  constants.  Hiki  and  Granato  used  pressure  derivatives  and 

uniaxial-stress  derivatives  to  determine  a  complete  set  of  TOE  constants 

g 

for  copper  at  room  temperature.  Salama  and  Alcrs  used  uniaxial  stress 
derivatives  exclusively  to  determine  a  complete  set  of  TOE  constants  for 
copper  at  three  different  temperatures,  295  °K,  77  °K,  and  4.2  °K. 
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Gaustcr  and  Brcazcalc  examined  combinations  of  copper  TOP  constants  at 
room  temperature.  Peters  and  Brcazcalc  and  Parc  extended  these 
measurements  to  77  °K.  In  this  paper,  we  report  results  of  copper 
which  have  been  measured  to  3  °K,  by  a  technique  sensitive  to  changes 
of  TOE  constants  as  a  function  of  temperature.  From  these  measurements, 
we  isolate  certain  combinations  of  TOE  constants,  which  are  of 
particular  theoretical  interest. 

The  noble  metals,  of  which  copper  is  an  example,  form  face-centered 
cubic  crystal  configurations,  for  which  simplified  models  exist  to 
explain  the  behavior  of  TOE  constants.  For  this  configuration,  we  find 
that  if  forces  of  interaction  arc  central  in  nature,  the  crystal  is  free 
from  external  stress,  and  each  atom  is  at  a  center  of  inversion,11  then 
the  Cauchy  relations  must  ho.d.' 


Second-order  constants 

C12  "  C44 

(1) 

Third-order  constants 

C1 12  =  C166 

(2) 

and 

f  =  r  =  r 
123  456  ^144 

g 

Hiki  and  Granato  have  shown  that  if,  in  addition  to  the  above 
assumptions,  nearest-neighbor  repulsive  interaction  is  the  predominant 
contribution  to  the  clastic  constants,  then  the  additional  relationships 
also  hold. 


Second-order  constants 

C11  " 

2C12  =  2C44 

(3) 

Third-order  constants 

Clll 

=  2Cn2  =  2C166 

(4) 

C1 23 

=  C456  =  C144  =  ° 

Our  data  allow  us  to  make  some  statements  about  the  validity  of 
the  TOE  constant  Hiki-Granato  relations  for  copper  between  3  °K  and  room 


temperature. 
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II.  EXPERIMENTAL  TECHNIQUE 


Pure  mode  propagation  for  a  longitudinal  ultrasonic  wave  is 

possible  for  three  principal  directions  in  a  cubic  crystal.  For  these 

12 

directions,  the  wave  equation  reduces  to 


(3K2  a  &  ^2 


(5) 


where  K2  and  are  combinations  of  SOE  and  TOE  constants  respectively, 
which  arc  given  in  Table  1. 

Assuming  a  sinusoidal  wave  of  frequency  0)  applied  at  a  =  0,  this 
equation  has  the  solution 

U  =  A^  sin(ka  -  cut )  -  (3K2  +  Kj)/8K0  A^2  k2a  cos  2(ka  -  cot;  +  . .  .  (6) 

where  k  is  the  propagation  constant  2tt/A,  a  is  the  propagation  distance 
in  the  sample,  and  A^  is  the  amplitude  of  the  fundamental.  The  amplitude 
A2  of  the  second  harmonic  term  is  given  by 

A2  =  -[(3K2  +  K3)/8K2]  Aj2  k2a  .  (7) 


The  measurement  of  A 2  and  Aj  is  the  basis  of  the  calculation  of  the 
nonlinearity  parameter  8  where  8  is  the  negative  of  the  ratio  of  the 
nonlinear  term  to  the  linear  term  in  Eq.  (5). 


8  =  - 


3K. 


(8) 


Solving  Eq.  (7)  gives  8  in  terms  of  measured  quantities: 


The  quantity  K.,  can  be  determined  by  the  relation  K2  =  pv*-,  where  v  is 
the  velocity  of  sound  in  the  appropriate  direction.  For  our  purposes, 
we  calculated  the  values  of  K2  at  the  various  temperatures  from  data 
given  in  Overton  and  Gaffney.^3 

III.  RESULTS  AND  DATA  ANALYSIS 

In  these  measurements  one  uses  techniques  similar  to  those  which  have 
been  previously  cited  in  the  literature. 3,^,3,6  The  room  temperature 
measurements  for  in  the  various  crystal lographic  directions  have 
been  taken  from  Peters,  Breazeale,  and  Pare,3  since  the  same 
samples  were  used.  Figure  1  shows  the  values  of  6  as  a  function  of 
temperature  in  the  three  principal  crystal lographic  directions.  Data 
from  Ref.  3  were  used  between  300  °K  and  77  °K.  Those  below  77  °K  arc 
the  new  data  which  were  matched  to  the  77  °K  datum  from  Ref.  3. 

Figure  2  shows  the  values  of  as  a  function  of  temperature 
calculated  from  the  data  of  Fig.  1.  The  scatter  in  the  [110]  data  results 
in  part  from  the  shape  of  the  sample.  The  [110]  faces  arc  at  an  angle 
of  approximately  15°  to  the  axis  of  the  cylindrical  sample.  This  canting 
of  the  sample  axis  led  to  difficulties  in  keeping  the  sample  seated  on 
the  ground  ring  of  the  capacitive  detector. 

Error  for  those  measurements  are  determined  by  the  measurements 
at  room  temperature  and  the  relative  measurements  at  the  other 
temperatures.  The  random  errors  for  at  room  temperature  arc  ±1.5", 

±3. 2°s,  and  ±20„  for  values  of  in  the  [100],  [110],  and  [111] 
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igure  1.  Measured  values  of  the  nonlinearity  parameter  as  a  function  of  temperature. 
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respectively.  Systematic  error  for  these  measurements  is  estimated  to 
be  at  most  ±10*.  It  is  estimated  that  K,  can  bo  measured  relative  to 
room  temperature  to  well  within  3%. 

Examination  of  Table  I  reveals  that  the  expressions  for  are  not 
the  simplest  combinations  of  TOE  constants  available  from  our  data.  The 
for  the  [100]  direction  is  the  single  TOE  constant  But  the 

for  the  other  directions  also  include  C^.  Thus,  it  is  possible  to 
subtract  from  for  the  directions  [110]  and  [111].  Proceeding 

in  this  fashion,  one  is  able  to  obtain  the  combinations 
and  C125  +  6C144  +  8C456  plotted  in  Fig.  3. 

IV.  DISCUSSION 

The  nonlinearity  parameter  B  in  Fig.  1  is  observed  to  be  relatively 
independent  of  temperature,  as  was  originally  assumed  to  be  the  case  with 
the  Griinclsen  parameter  y.1^'15  As  a  matter  of  fact  it  is  possible  to  make 
a  specialized  definition  of  a  "Gruneiscn  number"  which  is  related  to  the 
nonlinearity  parameter  ^  Variation  in  the  value  of  B  docs  occur  between 

approximately  25  °K  and  50  °K,  with  the  most  distinct  variation  occurring 
in  the  data  for  the  [110]  direction. 

When  one  examines  the  third-order  elastic  constants  combinations 
plotted  in  Fig.  2,  the  temperature  dependence  becomes  more  pronounced. 
Although  the  K3  for  [110]  direction  Cjjj  varies  by  only  8*  over  the 
temperature  range,  and  the  value  of  at  0  °K  is  almost  identical  to 

the  value  at  300  °K,  the  other  two  orientations  exhibit  more  variation 
with  temperature.  Nevertheless,  it  may  be  worthwhile  to  point  out  that 
the  three  curves  behave  in  somewhat  the  same  way.  The  most  obvious 
temperature  variation  occurs  in  the  Kj  for  the  [110]  direction.  The 


Figure  3.  Simplest  combinations  of  TOE  constants  plotted  as 
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origin  of  this  effect  can  he  located  somewhat  more  exactly  by  examining 
Fig.  5,  a  plot  of  the  simplest  TOC  constant  combination  available  from 
our  data. 

The  combinations  of  TOE  constants  in  Fig.  3  happen  to  correspond 
to  combinations  which  occur  in  the  Hiki-Granato  relations.  One  finds 
that  the  Hiki-Granato  relations  predict  that 

C 1 1 2  +  4C166  =  5/2  C1U 
and  C123  +  6C144  +  8C456  =  0. 

Thus,  on  Fig.  3  we  have  plotted  5/2  to  aid  in  the  comparison.  In 

Fig.  3  one  finds  that  Cjp  +  40^^  =  5/2  C  ^  to  within  approximately 
6"o  over  the  entire  temperature  range.  We  may  also  point  out  that  the  two 
curves  have  almost  identical  shapes  over  the  entire  temperature  range, 
and  this  implies  that 


3C, , ,  +  40,^)  3C 


'111 


'166 J 
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over  the  same  temperature  range.  This  observation  is  consistent  with 
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the  contentions  of  Hiki,  Thomas,  and  Granato  that  higher-order 
elastic  constants  of  materials  which  have  markedly  overlapped  closed 
shells  are  influenced  most  strongly  by  nearest  neighbors. 

The  combination  of  TOE  constants  +  6<“144  +  *^456^  exhibits 

an  interesting  behavior,  becoming  slightly  positive  above  200  °K. 

(This  combination  should  be  zero  according  to  the  central  force, 
nearest  neighbor  model.)  At  all  temperatures,  it  remains  small  in 
comparison  to  the  other  combinations.  However,  it  exhibits  a  dip  which 
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begins  near  45  °K,  reaches  its  minimum  value  near  75  °K  and  changes 
slope  in  the  neighborhood  of  200  °K.  Thi ,  behavior  has  the  same  general 

in 

temperature  dependence  as  the  Bordoni  peak  in  copper.-  Perhaps  this 
combination  is  sensitive  to  dislocation  movement. 

In  conclusion,  we  feel  that  our  data  are  nominally  consistent 
with  the  predictions  of  a  central  forces,  nearest  neighbor  interaction 
model.  As  usual,  there  are  details  which  need  to  be  explained,  but 
further  explanation  would  depend  upon  a  more  detailed  model  than  we 
have  used,  and  more  detailed  data  than  are  available. 
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We  present  results  of  examinat  :or  of  'he  u]  t  rasc-r.  ;c  harncru 
generation  technique  for  measurerent  of  nor  1. near  it  v  parameter*, 
and  heme  the  thir4  order  elastu  const  a:  t  *  .  of  solids.  Hereto¬ 
fore.  samples  cir  in  diameter  and  at  leas?  .  cm  in  length 

have  been  used  in  measurement  s  with  a  1  or.  diameter  V  HH:  ouart: 
transducer.  Our  results  reveal  that  the  mirmum  *>i:e  of  a  sample 
which  will  allow  accuracies  between  :  and  l''4  '  r,  .  :}y  measure¬ 
ments  is  4  mr  in  ultrasonic  path  length  and  ^  rr  in  diameter 
Correction  of  the  data  for  diffraction  effects  i«  necessary  with 
transducers  or  receivers  smaller  than  1  ~  .  *  5  •'jr. . 


INTkPIM-t’TIOV 


A  finite  amplitude  sinusoidal  ultr.'onic  wave  propagating  through  a  nonlinear  medium 
distort*  as  it  propagates.  This  waveform  distortion  t.»n  be  characterized  as  har¬ 
monic  generation,  and  the  measurement  of  the  harmonic  generation  makes  possible  the 
evaluation  of  the  nonlinearity  parameter*  of  the  material.  If  the  mater. al  is  a 
cubic  crystal,  then  one  finds  pure  node  propagation  along  the  principal  directions. 
Measurement  of  the  harmonic  generation  along  the  pure  mode  directions  reduces  to  the 
measurement  of  the  amplitude  of  *he  fundamental  A]  and  of  the  second  harmonic  As, 
since  the  nonlinearity  parameter  is  given  by  [1] 

f.«,  •  1.^1  SA, 


(  :  j  a;  r. 


where  and  Kj  are  'he  linear  combinations,  respectively,  of  the  second-  and  third 
order  e.astic  constants  given  in  Table  !.  k  =  is  the  propagation  constant,  and 
a  is  the  propagation  distance  sample  length'' 


Table  I  A,  and  k.  for  [1  (><'[,  [IK'1,  and  [1111  direction*  in  cubic  trs-tal* 
direction  >  ,  k  ^ 
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1  Xl'I  R TM*  NT  \l  IF  OiM;?n 

The  amplitudes  of  the  fundamental  and  of  the  second  harnonit  ■  measure.,  '■v  use  of 
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the  capacitive  receiver  shown  in  Fig.  1.  Ultrasonic  wave  pulses  are  generated  by  the 
quartz  transducer  and  propage  downward  through  the  sample  and  impinge  on  the  bottom 
surface.  Since  the  sample  surfaces  are  optically  flat,  a  receiver  button  can  be 
placed  at  a  distance  Sq  from  the  sample  to  form  a  parallel  ’late  capacitor.  Kith  a 
bias  voltage  of  approximately  ISO  volts  dc ,  one  finds  that  the  ultrasonic  wave 
impinging  on  the  interface  gives  an  ac  voltage  (.'] 


where  A  is  the  amplitude  of  the  harmonic  to  be  measured. 

Approx imat ions 

We  have  examined  the  approximations  inherent  in  the  use  of  harmonic  generation  to 
measure  the  nonlinearity  parameters  of  solids  and  find  that  three  of  then  need  to  be 
examined  in  order  to  define  the  smallest  sample  one  can  use 
(a^  Infinite  plane  wave  assumption, 

(M  Parallel  plate  capacitor  approximation, 
fc'l  Iffect  of  diffraction. 


ASAl>c  IS  Of  MI  A^UM  Ml  NT  T!  OINIQUf 

With  the  approximations  in  mm.,  one  can  evaluate  the  smallest  sample  ^ize  us., Me 
with  present  techniques- 


a  '  Infinite  P 1 have  Assumption 


In  deriving  5q.  il),  the  assumption  is  mac.i  that  the  record  harmonic  i c  generated  by 
a  fundamental  which  is  an  infinitely  extended  plane  wave.  Although  thi*  m  not 
strictly  true  in  our  experimental  arrangement,  as  we  will  see,  we  can  assume 
infinitely  extended  plane  waves  for  the  moment  to  estimate  the  minimum,  sample  ?  h ; c  K 
ness  required  to  produce  measurable  second  harmonics.  We  have  found  that  we  car 
reliav!>  measure  second  harmonics  Aj*  of  the  order  of  Angstrom* ,  and  have 

measured  lf  1  Angstroms.  Using  A  "*ir  -  !<•'  **  Angstrom*,  and  value*  of  I;  and  »  ;  frr 
copper  [11!  1  as  t  epical  values,  nr.**  finds  that  the  mrmc  sample  thicknc**  arir 
d  mm.  In  ■'♦her  material*  different  value*  of  *  ■  and  *  -s  ".it?  charge  th;*  -irinr 
thi.iress  b\  a  f  ictor  ot  two  ir.  either  direction 
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was  assumed  to  satisfy  the  infinite  plane  wave  assumption  reasonably  well  in  samples 
which  were  small  enough  to  place  in  a  helium  cryostat. 


EXPFRIMENTAL  RESULTS 

We  now  would  like  to  report  the  results  of  a  systematic  experimental  investigation  of 
the  effect  of  diffraction  on  our  measurements.  A  single  crystal  copper  sample  was 
selected  which  allowed  ultrasonic  wave  propagation  along  the  [111]  direction.  The 
sample  was  2.54  cm  in  diameter  and  3.96  cm  in  length.  By  selecting  transducers  with 
diameters  given  in  Table  II,  we  were  able  to  effectively  have  the  capacitive  receiver 
in  either  the  Fresnel  zone  or  the  Fraunhofer  zone  of  the  ultrasonic  wave  diffraction 
field.  Two  capacitive  receivers  were  used  with  diameters  of  11.8  and  6.36  mm. 

Table  II.  Transducer  diameters  and  ultrasonic  beam  characteristics  used  in  a  [111] 
copper  sample.  Receiver  diameters  were  11.8  and  6.36  mm. 


Transducer 

Diameter 

(mm) 

d/> 

a*7> 

(mm) 

g 

(degrees) 

12.65 

74 

233 

0.9 

5.40 

32 

43 

T  n 

5.’0 

:o 

3.2 

1 . 95 

11 

6 

6.1 

To  give  an  idea  of  the  relative  size,  a  scale  drawing  of  the  (linear]  ultrasonic 
wave  field  in  the  sample  is  given  in  Tig.  2.  To  make  the  drawing  we  have  assumed 
that  the  intermediate  size  3 . 20  mm  transducer  is  behaving  as  a  piston  vibrator  and 
is  producing  a  beam  whose  half-angle  0  is  given  by 

sin  B  *  0.(1  )/a  (3) 

where  >  is  the  (fundamental)  ultrasonic  wavelength  and  a  is  the  transducer  radius. 

The  distance  a*7\,  sometimes  referred,  to  as  the  Tresne1.  distance,  also  is  indicated. 
For  the  largest  transducer  this  distance  is  much  greater  than  the  sample  length  so 
that  the  measurements  are  made  within  the  Tresne!  zone  where  the  plane  wave  approxi¬ 
mation  l*  reasonably  we  1 1  - sa t i sf i ed . 

As  seen  from  Fq.  (I),  a  measure  of  the  nonlinearity  parameter  for  a  given  frequency 
and  sample  length  is  the  ratio  From  Fq.  (2',  it  is  clear  that  this  ratio  is 

proportional  to  \  VVj  where  V;  and  are  the  voltages  generated  in  the  receiver  by 
the  fundamental  and  the  second  harmonic,  respectively.  We  would  determine  the  effect 
on  this  measured  a*antitv  b>  diffraction.  The  averages  of  a  number  of  measurements 
of  the  ratio  V'Ar  are  plotted  as  a  function  of  the  Transducer  diameter  in  big.  3. 

As  indicated,  we  find  that  the  -e.i*..rod  '  j  A  j  values  fall  or,  smooth  curves  for  each 
re  eiver  size.  At  the  moment  there  i*  no  theoretical  justification  for  this  experi¬ 
mental  lv  observer  fact,  however,  '  ere  are  some  „  onM'tent  observations  which  can 
Ke  made. 

The  :at:  r-'irt  f  t  the  t  r  jro  ii.«.  er  •  r  ece  l »  er  i  '»"-b  i  r  » t  :  or  heretofore  used  in  mcu-'iir  i  ng 
rv  «r  I  i  r.o.i  r  1 1  •'  nrv*e*crr  :  s  indicate. I  by  .♦  so  i  i  f  point.  It  is  re.i'-urirg  to  note  that 

V*  is.-  j  rr'tc  -f  • v  e  v  .  r.e  f  sr  ■  ;  e  "iller  re  ''i.er  passes  Through  tin-  <>J  ul 

r  •.-•  i;m  ac  f  the  i  •■ter,  .  •'  *'-n  lata,  ":c‘*  vru~e  that  out  .bit  a  to 

*»'e  ‘  a  .  *•  •  t  i*  er.  .rj'-  .  c.r.d:ti  "  -  ^  :.1'  rr.vrr  il  I,  wri;  sat  i  s  f .  t  •  »•  t  ret  1 1  a  1 

*  .  i*  :  thus  •->  r  t  rea  ; .  r  c  .  -  rre.  t  i  on  .it  ? !: :  *  ;**  e  .  I  art  bet  ,  the  con- 

;  •!(■'  .  •  f'e  ;  r  •.  e «  v  .  g  c  •*  *  *  *  1  *  c~u.lcr  a”-: 1  <"  v.iT.  ■  e  .;*<•!  than  f>  have  used 

•  frt’f.  *  I  r  ..it  a  .  v  1  **  *  in  -a-p  I  <-s  s.-  s»  a  ;  ;  •  *  :t  the  plant  wave 

i*  i-  -at  Ir  '  a  i  • .  c  ur\  e<  m  - :  1 ..  r  *  '  ive  m  Mg.  3  car. 

.  -  •  •«  i*  riti  •?  .  r*  to  i '. . .  1 1  e  •re-  !  a  t  a  taker,  wi*  ~  a  .  1  sample-  the 

■■ :  ;  v-d  her.  e  t  r  - o-  .  ;  era r  ;  •  •  r  ir.v*eter«  '  A  ~  'ic  ter  l  st  i ,  i 

M:,e  m  •  .» *  *.•  viTi'fv  t  ‘  ^  *' *' •  ~  * -*•  „ ;  me  wave  nr  r;;T‘  For  ( 


lo,  the 
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data  points  at  a  transducer  diameter  of  3.7  mm  require  a  correction  to  bring  them  in 
line  with  the  value  of  Vi/V^  indicated  by  the  solid  point,  and  the  correction  factor 
can  be  evaluated  from  the  curves.  Correction  of  data  for  smaller  transducer 
diameters  might  not  be  accurate  enough.  Allowing  the  sample  walls  to  be  well  out¬ 
side  the  ultrasonic  beam  might  require  a  sample  diameter  of,  say,  S  mm.  Thus,  we 
come  to  an  approximate  answer  to  the  question  po^cJ:  the  smallest  sample  usable  with 
the  harmonic  generation  technique  has  a  diameter  of  5  mm  and  a  length  of  4  mm.  The 
surfaces  must  he  optically  flat  and  parallel,  of  course. 

Finally,  we  would  point  out  that  a  very  promising  prospect  exists  for  progress  in  the 
general  subject  of  nonlinear  acoustics.  Figure  3  is  the  first  numerical  evidence  we 
have  that  after  including  the  nonlinear  terms  describing  the  propagating  medium,  the 
diffraction  integral  should  have  solutions  which  are  tractable.  The  curves  give  the 
values  of  an  integral  across  a  plane  at  specified  distances  from  a  sinusoidally 
vibrating  piston  source  which  is  radiating  into  a  nonlinear  medium.  A  large  number 
of  such  curves  taken  for  different  distances  and  different  transducer  sizes  would 
allow  one  to  evaluate  directly  the  effect  of  diffraction  on  harmonic  generation. 

This  remains  one  of  the  unsolved  fundamental  problems  of  nonlinear  acoustics. 
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The  relationship  between  the  ultrasonic  nonlinearity  parameter  for  solids  and  the  acoustic  Griineisen 
number  has  been  derived  for  longitudinal  ultrasonic  wave  propagation  in  the  pure  mode  directions  of 
cubic  crystals  and  isotropic  solids.  Agreement  between  the  acoustic  Griineisen  number  and  the 
thermodynamic  Griineisen  parameter  is  best  for  ultrasonic  harmonic  generation  measurements  along  the 
[100|  direction  of  a  cubic  crystal.  Comparison  of  the  temperature  curves  of  the  acoustic  Gruneisen 
number  of  copper  shows  that  the  acoustic  Gruneisen  number  generally  follows  the  temperature 
dependence  of  the  lattice  contribution  to  the  thermodynamic  Gruneisen  parameter. 

PACS  numbers:  43.25  Ba.  43.35.Cg 


INTRODUCTION 

In  developing  the  harmonic  generation  technique  for 
measurement  of  the  nonlinear  properties  of  solids,  we 
have  found  a  fundamental  significance  to  the  ratio  of 
coefficients  of  the  nonlinear  terms  to  the  linear  terms 
in  the  nonlinear  wave  equation  describing  the  propaga¬ 
tion  of  a  finite  amplitude  ultrasonic  wave  in  the  solid. 
This  ratio,  defined  as  the  ultrasonic  nonlinearity  pa¬ 
rameter  for  solids,  is  found1  in  the  perturbation  solu¬ 
tion  for  the  particle  displacement,  in  the  expression 
for  the  discontinuity  distance,  in  the  implicit  solution 
for  the  wave  velocity,  etc.  In  addition,  this  quantity  is 
observed  to  be  only  weakly  dependent  on  temperature 
in  those  solids  studied. 2-5  These  observations  have  led 
us  to  realize  that  there  is  an  intimate  relationship  be¬ 
tween  the  ultrasonic  nonlinearity  parameter  and  the 
Gruneisen  parameter  evaluated,  for  example,  from 
studies  of  Brillouin  scattering,8  from  thermal  expan¬ 
sion,7  or  from  ultrasonic  attenuation.8  The  purpose  of 
the  present  work  is  to  define  the  acoustic  Gruneisen 
number  for  cubic  crystals  and  isotropic  solids  and  to 
give  the  relationship  between  it  and  the  ultrasonic  non¬ 
linearity  parameter.  Values  of  the  acoustic  Gruneisen 
number  taken  from  room  temperature  data  on  ultra¬ 
sonic  harmonic  generation  are  presented,  and  a  com¬ 
parison  is  made  between  the  temperature  dependence 
of  the  acoustic  Gruneisen  number  and  that  of  the  ther¬ 
modynamic  Gruneisen  parameter. 


I.  THEORY 

One  can  consider  the  wave  equation  in  a  cubic  crys¬ 
tal  for  three  pure-mode  directions:  ( 100 1,  (110 1,  and 
(111).  In  these  three  directions  pure  longitudinal 
waves  may  propagate  and  the  transverse  wave  is  not 
excited.  The  nonlinear  equation  of  motion  for  pure  lon¬ 
gitudinal  waves  in  these  three  directions,  assuming  no 
attenuation,  may  be  written 
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where  p0  is  the  unstrained  mass  density,  tt  is  the  lon¬ 
gitudinal  displacement,  and  a  is  the  Lagrangian  coor¬ 
dinate  along  the  wave  propagation  direction.  A',  and  A’3 
are  linear  combinations  of  second-  and  third-order 
elastic  constants  and  are  given  in  Table  I  for  the  three 
pure-mode  directions  in  the  cubic  crystal.  The  equa¬ 
tion  of  motion  for  an  isotropic  solid  is  identical  to  that 
of  a  cubic  crystal  in  the  (100]  direction,  with  the  ap¬ 
propriate  interpretation  of  Cu  and  Cm.  An  implicit 
solution  for  the  particle  velocity  3tr/3/,  which  satisfies 
Eq.  (1)  and  the  boundary  condition 


57-(5L“""'  ® 

can  be  written  in  the  form1 


3 n  (hu\  (  ,  f/A'A 

—  =l — I  sinttW  -  tea  (  — 

3'  \a/  /„«,  \  LW 

l/A'A"1 /3A,  +  A'3\ 3nl  I 

+  2\Po )  V  l<2  'Mr 


(3) 


In  analogy  with  the  results  of  Earnshaw9  for  gases,  one 
finds  that  the  phase  velocity  Vf  can  be  written 


vr-v0i*\^t 


(4) 


where  V0(=  (A2/p0)1  /2  is  the  velocity  of  a  wave  having 
infinitesimal  amplitude  and 


TABLE  I. 

A’;  and  A's  for  IlOO), 

1110),  and  llll)  directions. 

Direction 

a-2 

a3 

liool 

C||l 

moi 

Cu  1  C  12  4  2  Cm 

Cm  4  3Cn2  4  12Ci66 

2 

4 

lull 

t  il  +  2Ci2  4  4 Cm 

Cm  4  6 Cm  <  lSCm  4  24t  i£C 

3 

9 

2C  123  4  16C4SC 
9 
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TABLE  II.  Acoustic  Gruneisen  numbers  V*  for  [100),  ( 110].  and  illl)  directions  in  cubic  crystals. 


Direction 


[1001 
1 110! 
Hill 


yf 


r  ' 

3Cn+Cui 

f  12C44  +  Cm  +  3Cii2  + 

2CU  J 

f 

6Cn  +6C12 

l  4C11+4C1j  +  8C44  J 

(9C11  +  I8C12  +  36C44  +C'm  +  6C112  +  12CU4  +24Ci66  +2C123 +  16C'45sl 

l 

6Clt  +  I2Cw  +  24C 

«  J 

+*,)/*„],  (5) 

is  defined  as  the  ultrasonic  nonlinearity  parameter  for 
solids.  The  subscript  i  represents  all  indices  specify¬ 
ing  longitudinal  wave  propagation  in  the  pure-mode  di¬ 
rections  of  the  crystal.  The  effect  of  the  nonlinear 
terms  in  the  wave  equation,  then,  is  to  change  the 
phase  velocity  by  an  amount  proportional  to  the  product 
of  the  particle  velocity  du/dl  and  the  nonlinearity  pa¬ 
rameter  &{.  This  solution  is  valid  for  propagation  dis¬ 
tances  less  than  the  discontinuity  distance 


which  is  of  the  order  of  1  m  for  common  solids.  From 
Eq.  (4)  and  assuming  irrotationality,  the  change  in  ve¬ 
locity,  AFj  =  resulting  from  the  nonlinearity 

of  the  medium  can  be  expressed  in  terms  of  the  dis¬ 
placement  gradient  3«/3 a  by 


1 

2 


3  u 

oi  3a  ' 


(7) 


In  order  to  relate  this  change  in  velocity  to  the 
change  of  thermal  phonon  velocity  resulting  from  lat¬ 
tice  anharmonicity,  it  is  convenient  to  define  the  acous¬ 
tic  Gruneisen  number  specifying  the  adiabatic  strain 
dependence  of  a  lattice  frequency  v,  for  the  longitudinal 
mode  i  by 


where  the  subscript  S  emphasizes  the  fact  that  the  de¬ 
rivative  is  taken  under  isentropic  conditions,  and  the 
subscript  i  represents  all  indices  specified  by  longitu¬ 
dinal  wave  propagation  in  the  directions  [100 J,  [110], 
and  [ill  j.  Here  we  take  the  derivative  with  respect  to 
the  displacement  gradient  3n/3 a  rather  than  with  re¬ 
spect  to  the  Lagrangian  strain  measure  as  ordinarily 
is  done  in  defining  generalized  Gruneisen  parameters. 


This  is  possible  because  in  harmonic  generation  exper¬ 
iments  the  nonlinearity  parameter  is  obtained  by  ex¬ 
trapolating  to  zero  amplitude  a  plot  of  the  ratio  of  the 
amplitudes  of  the  second  harmonic  to  the  square  of  the 
fundamental  amplitude  as  a  function  of  the  fundamental 
amplitude.  In  the  limit  of  zero  amplitude  the  difference 
between  the  displacement  gradient  and  the  Lagrangian 
strain  measure  vanishes.  When  evaluated  at  zero 
strain,  the  acoustic  Gruneisen  number  of  Eq.  (8)  is  the 
same  as  Brugger’s  tensorial  isentropic  Gruneisen  num¬ 
ber10  for  a  longitudinal  strain  along  the  propagation  di¬ 
rection  of  the  mode. 

According  to  the  standing  wave  condition  of  the  Debye 
continuum  model,  for  any  state  of  strain  the  mode  fre¬ 
quencies  are  related  to  the  wave  speed  V,  and  the  un¬ 
strained  dimension  l  of  the  crystal  by10 

vtoc  Vjl.  (9) 

Substituting  this  expression  for  the  mode  frequencies 
into  Eq.  (8)  and  integrating  between  Voi  and  Vjt  one  ob¬ 
tains 

(10) 

Comparing  Eqs.  (10)  and  (7)  one  finds  the  relationship 
between  the  acoustic  Gruneisen  number  and  the  ultra¬ 
sonic  nonlinearity  parameter  for  solids  is  expressed 
by 

rf  =  i0,.  (11) 

In  Table  II  are  listed  (he  acoustic  Gruneisen  numbers 
for  the  three  pure-mode  directions  in  cubic  crystals. 

II.  COMPARISON  WITH  EXPERIMENT  AND 
CONCLUSION 

Ultrasonic  nonlinearity  parameters  have  been  mea¬ 
sured  for  germanium,  copper,  and  fused  silica.  As 


TABLE  III.  Comparison  of  adiabatic  Gruneisen  number  and  other  Gruneisen  parameters. 


yf  r 

Present  work  Nava  and  Romero-1 

rfiool  rfnoi  tfuil  1  f  loo]  l|uo|  1  (ml 


Copper  2.64  5.50  3.94  2.00 

Germanium  1.5  3.0  2.6  1.07  1.12  1.41  0.75 

Fused  silica  —5.x  O.ix 


V 

Collins  and  White1' 


'  It.  Nava  and  J.  Romero,  J.  Aeoust.  Soc.  Am.  64.  529—532  (1978). 

h  J.  G.  Collins  and  G.  K.  White,  “Thermal  Expansion  of  Solids"  in  Progress  in  I. on  Temperature 
1'lnisics,  edited  by  C.  J.  Gorter  (Wiley.  New  York.  1964),  Vol.  4. 
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these  measurements  have  been  made  from  room  tem¬ 
perature  down  to  3  °K,  a  fairly  detailed  comparison  can 
be  made  between  the  acoustic  Griineisen  number  and 
other  Griineisen  parameters  in  some  cases. 

In  Table  III  we  list  room  temperature  values  of  the 
acoustic  Griineisen  number.  Also  listed  are  room  tem¬ 
perature  values  of  thermodynamic  Griineisen  param¬ 
eters’  and  effective  ultrasonic  Griineisen  parameters. 
The  thermodynamic  Griineisen  parameter  is  given  by 

y=a/K1Cv=a/KsCp,  (12) 

where  a  is  the  total  thermal  volume  expansivity,  KT 
and  Ks  are  the  isothermal  and  the  isentropic  compres¬ 
sibilities,  respectively,  and  CY  and  C p  are  the  isocho- 
ric  and  isobaric  heat  capacities.  The  effective  ultra¬ 
sonic  Griineisen  parameter  Fp  is  related  to  the  ultra¬ 
sonic  attenuation  for  a  wave  of  polarization  p  and  prop¬ 
agation  direction  q  by” 

af=(3n2A',r/Pc#s2)r>,  (13) 

where  A'a  is  the  thermal  conductivity  along  q,  T  is  the 
absolute  temperature,  Q  is  the  ultrasonic  angular  fre¬ 
quency,  p  is  the  mass  density,  and  C „  and  S  are  the 
sound  wave  and  Debye  average  velocities,  respectively. 

Examination  of  Table  III  reveals  the  fact  that  the 
acoustic  Griineisen  number  yf,no)  agrees  most  closely 
with  the  thermodynamic  Griineisen  parameter,  and  that 
all  three  of  the  acoustic  Griineisen  numbers  are  larger 
in  magnitude  than  the  thermodynamic  Griineisen  pa¬ 
rameter.  The  greatest  discrepancy  between  the  two 
Griineisen  parameters  occurs  for  fused  silica  for 
which  the  acoustic  Griineisen  number  is  negative.  In¬ 
terestingly  enough,  the  behavior  is  consistent,  for 
fused  silica  exhibits  an  anomalous  behavior  of  many 
of  its  thermodynamic  and  ultrasonic  properties. 12'14 
The  negative  value  of  the  acoustic  Griineisen  number 
is  related  to  the  fact  that  the  second  harmonic  in  fused 
silica  is  generated  180'"  out  of  phase  with  that  generated 
in  other  solids.15 

It  is  significant  that  the  agreement  between  the  acous¬ 
tic  Griineisen  numbers  and  the  effective  ultrasonic 
Griineisen  parameters  of  germanium  is  as  good  as  it  is, 
because  the  data  come  from  fundamentally  different 
measurements:  ours  from  harmonic  generation,  Nava 
and  Romero’s  from  attenuation. 

The  temperature  dependence  of  the  acoustic  Griinei- 
sen  number  for  the  [100]  direction  in  single  crystal 
copper  was  found  to  be  quite  similar  to  that  of  the  ther¬ 
modynamic  Griineisen  parameter.  The  two  are  plotted 
on  the  same  graph  in  Fig.  1.  The  temperature  scale  is 
the  reduced  temperature  T/8,  where  the  Debye  tem¬ 
perature  8  =  345  K  for  copper.  A  smooth  curve  is 
drawn  through  experimental  values  of  y5llKI].  For  com¬ 
parison,  values  of  the  thermodynamic  Griineisen  pa¬ 
rameter7  are  plotted  on  the  same  graph.  The  total 
Griineisen  parameter  is  indicated  by  a  solid  curve.  In 
the  low  temperature  region  the  electronic  contribution 
to  the  thermodynamic  Griineisen  parameter  is  distin¬ 
guishable.  The  dotted  curve  labeled  ‘7"  for  “lattice" 
contribution  remains  essentially  constant,  while  the 
total  carve  labeled  ‘7  +  c"  for  “lattice  plus  electron" 
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FIG.  1 .  Temperature  dependence  of  the  acoustic  Griineisen 
number  yf,ool  and  the  thermodynamic  Griineisen  parameter  y 
of  copper. 

decreases  as  the  temperature  is  lowered. 

It  is  interesting  to  notice  that  the  acoustic  Griineisen 
number  yf,001  remains  parallel  to  the  lattice  contribu¬ 
tion  to  the  thermodynamic  Griineisen  parameter,  and 
almost  exactly  1.5  times  as  large.  The  significance 
of  this  observation  is  under  investigation  at  the  present 
time. 
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ABSTRACT 

Nonlinear  interaction  of  ultrasonic  waves  is  discussed  from  the 
very  general  quantum  mechanical  viewpoint.  By  using  the  concept  of  three- 
phonon  interactions  one  is  able  to  derive  formulae  for  power  and  intensity 
of  ultrasonic  waves  generated  through  nonlinear  mixing  of  two  ultrasonic 
waves,  as  well  as  the  formula  for  parametric  amplification  of  an  ultra¬ 
sonic  wave  through  its  nonlinear  interaction  with  another  ultrasonic  wave. 
The  effect  of  attenuation  on  nonlinear  ultrasonic  wave  interactions  also 
is  discussed,  and  it  is  shown  that  in  some  situations  attenuation  strongly 
affects  the  generation  of  the  second  harmonic  of  an  initially  sinusoidal 
ultrasonic  wave.  It  is  shown  that  in  the  correspondence  limit  the  quantum 
mechanical  viewpoint  gives  results  which  arc  in  very  good  agreement  with 
results  obtained  from  classical  physics. 


Will  be  published  in  November  1980  issue  of  Journal  of  the  Acoustical 
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I.  INTRODUCTION 

Nonlinear  interactions  of  ultrasonic  waves  in  solids  arc  of  two  general 
types.  First  is  the  interaction  in  which  two  waves  generate  a  third  wave 
(this  also  includes  harmonic  generation  of  an  initially  sinusoidal  wave,  as 
we  will  sec).  Second  is  the  amplification  of  an  ultrasonic  wave  through 
its  interaction  with  another  ultrasonic  wave.  Taken  together,  these  two 
types  of  nonlinear  interaction  offer  a  wide  range  of  prospects  of  technical 
application  in  addition  to  providing  an  excellent  possibility  to  study  the 
fundamental  anharmonicity  of  a  crystalline  lattice. 

Theoretical  discussion  of  the  subject  can  be  either  from  the  classical 

approacn  or  the  quantum  mechanical  approach.  In  the  classical  approach  one 

uses  coupled-mode  equations*  6  and  the  classical  language  of  "plane  waves 

propagating  in  a  semi-infinite  medium."  In  spite  of  the  fact  that 

classical  physics  is  quite  adequate  for  discussion  of  ultrasonic  interactions 

in  which  quantum  effects  are  negligible,  the  quantum  mechanical  approach 

7-12 

using  the  concept  of  phonons  also  has  been  applied  to  this  problem.  The 

quantum  mechanical  approach  gives  a  reasonably  straightforward  solution  even 
when  the  interaction  is  treated  as  a  three-dimensional  problem  in  an 
anisotropic  medium. 

In  general,  the  advantages  of  the  quantum  mechanical  point  of  view  in 
the  description  of  ultrasonic  wave  interactions  have  not  been  utilised  fully. 
Occasionally  discrepancies  occur  in  the  results  obtained  by  different 
authors,7’**  and  often  there  is  difficulty  in  comparing  different  theoretical 
results  because  the  definitions  of  the  physical  parameters  arc  not  totally 
consistent.  Such  nonuniformity  of  the  theoretical  approach  tends  to 
discourage  comparison  between  experimental  results  on  the  nonlinear  behavior 
of  ultrasonic  waves  and  results  of  quantum  mechanical  theories. 


i 


The  present  paper,  therefore,  is  devoted  to  a  unified  quantum 
mechanical  discussion  of  the  nonlinear  behavior  of  ultrasonic  waves  in 
crystals.  We  discuss  both  noncollincar  and  cnllineur  interactions  of  ultra¬ 
sonic  waves  and  point  out  the  possibility  to  measure  nonlinear  constants  of 
the  propagating  medium.  We  employ  the  quantum  mechanical  treatment  using 
the  concept  of  phonons  and  show  that  the  quantum  mechanical  treatment  gives 
not  only  results  obtained  previously  by  classical  methods,  but  also  gives 
the  possibility  to  describe  very  general  cases  of  interaction  of  ultrasonic 
waves  in  an  anisotropic  medium  of  finite  size. 

In  describing  the  nonlinear  behavior  of  ultrasonic  waves  using  the 
concept  of  three-phonon  interaction,  the  very  first  problem  is  the  derivation 
of  the  interaction  Hamiltonian  density  operator.  Section  II  is  devoted  to 
this  derivation.  Section  III  dcscriues  interactions  in  which  two  ultrasonic 
waves,  through  nonlinear  mixing,  generate  an  ultrasonic  wave  whose  frequency 
is  the  sum  of  the  frequencies  of  the  mixed  waves  (parametric  up-conversion). 
Section  IV  discusses  parametric  down-conversion  of  two  ultrasonic  waves 
and  parametric  amplification  of  an  ultrasonic  wave  through  its  interaction 
with  another  ultrasonic  wave.  In  Section  V  we  consider  the  effect  of 
attenuation  in  nonlinear  interactions  of  ultrasonic  waves.  Finally,  in 
Section  VI  we  discuss  the  effect  of  our  limiting  assumptions  and  relate 
them  to  assumptions  made  in  classical  theories. 
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II.  INTERACTION  MAM  f  l/l'ON  I  AN 


The  interaction  Hamiltonian  density  operator  for  three-phonon  interactions 

can  be  obtained  by  starting  with  the  classical  form  of  the  energy  of  interaction 

between  two  ultrasonic  waves.  First,  one  replaces  the  classical  displacement 

vector  by  the  corresponding  operator.  The  components  of  the  displacement 

1  1 

operator  then  have  the  form: 


u. 

i 


2pVw  J  1  q 


iq*r  „ 
e  +  a* 


-iq*r 


(q) 


(i) 


where'll  is  Planck's  constant  divided  by  2tt;  p  is  the  mass  density;  V  is  the 
volume  of  normalization;  wq  is  the  angular  frequency  of  a  phonon  having  the 
wave  vector  q;  a  and  a*  are  respectively  the  annihilation  and  the  creation 

q  q 

operators  of  the  considered  phonons;  and  are  components  of  the  polarization 

vector  (q) 

k™  -  ^ 


(q)  ’ 


Cm  1 

where  UQ  1  is  the  amplitude  of  the  ultrasonic  wave.  As  the  phonon  of  wave 
vector  q  in  general  can  have  three  different  polarizations,  the  index  "q"  will 
be  understood  to  be  a  double  index,  referring  both  to  wave  vector  and  to 
polarization.  Here,  and  in  what  follows,  the  Einsteinian  summation  over 
repeated  indices  will  be  understood. 

Next,  the  classical  form  of  an  interaction  energy  density  for  the  scattering 

of  an  ultrasonic  wave  S. .  by  an  ultrasonic  wave  Sf .  can  be  obtained  if  we  assume 

ij  ij 

that  the  wave  Sf.  modulates  the  properties  of  the  medium  in  which  the  wave  S .  . 
is  propagating.  Thus,  the  energy  density  of  the  wave  has  the  form 
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U  = 


^ iCi jkfi.  + 


3C.  ..  „ 
ljkJl 

3S 

mn 


r  ]  s. . j 

mn  J  1  j 


ki 


(2) 


where  the  S. .  and  the  Sf.  are  the  strain  tensor  components  of  the  given 
ij  lj 

ultrasonic  waves  and  the  are  t^e  elastic  moduli  of  the  medium  stiffened 

by  the  presence  of  an  ul trasonic  wave. 

Using  relations  (1)  and  (2)  and  making  a  plane  wave  expansion,  brings  the 
Hamiltonian  density  to  the  form 


h  - 

int  ~  ~  2  (2pVj 


ijka  \  j  k  l 


fk 


*(q) 
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q  q  J{  q  q 


^  _  '-t 
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V' 


where  v  ,  v ^  and  v^>  are  phase  velocities  of  the  phonons  (q,o^),  (q',u)^„)  and 

(q" respectively,  and  rafq^ ,  mjq  ^ ,  and  mjq  ^  are  components  of  the  unit 

vectors  n/q^  =  q/q,  m^q  ^  =  q'/q'  and  m^q  ^  =  q"/q".  The  first  part  of  the 

interaction  Hamiltonian  density  (3)  arises  as  a  result  of  the  nonlinearity  of 

the  strain  tensor  and  it  expresses  self-coupling  of  phonons;  the  second  part  of 

14 

(3)  is  a  consequence  of  the  anharmonicity  of  the  lattice.  In  this  form,  the 


Hamiltonian  density  describes  two  kinds  of  three-phonon  processes:  First,  the 
creation  of  a  phonon  through  the  annihilation  of  two  phonons  is  suggested  by 
those  terms  in  Eq.  (3)  which  contain  a^,  8^,,  and  8*,„;  second,  the  disintegration 


4b 


of  a  phonon  and  the  siinultaneou'  creation  of  two  new  phonons  is  suggest ed  by 
terms  containing  .1  ,  3*^,  ami  3*>... 

q  q  *  q 

In  the  following  we  discuss  both  processes  and  show  their  macroscopic 
behavior  in  the  nonlinear  interactions  of  ultrasonic  waves. 

III.  PHONON  CRLATION  THROUGH  ANNIHILATION  OF  TWO  PHONONS 
(Nonlinear  Parametric  Up-Conversion) 

In  a  nonlinear  medium  two  phonons  can  interact  to  create  a  new  phonon.  The 
experimental  situation  corresponding  to  this  process  is  the  situation  in  which 
two  ultrasonic  waves  of  angular  frequency  and  generate  a  new  ultrasonic 
wave  of  angular  frequency 

In  describing  the  main  features  of  this  interaction,  we  will  need  knowledge 
of  the  size  and  shape  of  the  interaction  volume.  In  order  to  keep  our  discussion 
general,  but  also  fairly  simple,  we  shall  assume  that  the  interaction  volume  and 

the  volume  of  normalization  coincide  with  the  volume  of  the  sample  under 

investigation.  The  interaction  coupling  constant  G  ,  however,  will  be 

*  1  6  qlq2q3 

assumed  to  have  nonzero  value  only  in  that  part  of  the  sample  in  which  the 

interaction  actually  takes  place.  This  assumption  can  be  employed  if  the 
considered  part  of  the  sample  is  much  larger  than  the  wavelength  of  the  ultra¬ 
sonic  wave. 

Let  us  now  assume  that  there  are  phonons  having  the  wave  vector  q^ 
and  N,  phonons  having  the  wave  vector  q,,  in  the  volume  of  the  sample.  Let  us 
also  assume  that  the  initial  state  of  the  newly  created  phonons  having  the 
wave  vector  q^  is  empty.  This  means  that  the  initial  state  vector  is 
|n  ,N7,0>.  The  probability  amplitude  for  transition  from  this  initial  state 
to  the  final  state  <Nj"l  * *  I  *s 
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wher 0  V  ,  the  interaction  volume,  is  the  volume  in  which  C  is  assumed 
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to  have  nonzero  value.  The  quantity  C.  ^  ^  ^  has  the  form 

f  (qJ  (q,)  (q3)  (q,)  _(q, )  Jq^) 
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where  (k  «k  “  )  is  the  cosine  of  the  angle  between  the  polarization  vectors 
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k  and  k  ,  etc. 

In  deriving  the  probability  amplitude,  Eq.  (1),  we  have  assumed  that 
N^>>1,  and  N^>>1.  A  factor  3!  is  implicit  in  Eq.  (4)  as  a  consequence  of  the 
summation  over  all  possible  q,  q' ,  and  q".  The  vectors  q^  and  q7  are  fixed 
by  the  experimental  situation,  and  the  transition  probability  amplitude  is  to 
be  considered  a  function  of  q  .  The  variability  of  q~,  however,  is  limited  by 
the  integral 
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The  function  .  depends  on  the  actual  sh..  <•  aia  >i:e  ,  t  *  he  ..  hunt-  ...  ♦ ! 

part  of  the  sample  m  which  interaction  actual lv  takes  place.  Therefore,  for 
a  detailed  treat  nent ,  the  value  of  '  should  l,t  calculated  foi  eas  h  e>p.er ’.rent 
situation  separately.  In  the  general  case  sue?-  et  t  >u  t  -  as  cor.u  .  1  re!  r;n  t  i .  >r. 
would  complicate  the  calculation  somewhat.  for  rare  ex;  er  ir.<  i.t  a  1  situations, 
however ,  it  is  adequate  to  proceed  as  tallows. 
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value  wile;,  ‘  q  /  l.i.,  for  q.  ♦  q,  i  q_. 

The  partial  transition  proh.il' i  !  1 1  .  rate  ot  the  t  ran -it  ion  un«.  i 
considerat  ion  can  now  lie  expressed  in  tin  form 
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♦  q ,  in,  respectively,  the  (x,yj  plane 
through  the  integral  V ,  which  is  expressed 
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=  e,  we  specify  that  |q^|  is  a  constant 

value  is  determined  by  sirr.ultaneous  fulfillment  of  the  law  of  conservation  of 

am!  the  lu»  >f  um  ervat ior  of  momentum  of  the  interacting  quasi - 
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In  tin-  following  we  .hall  iliscn  some  special  examples  of  ultrasonic 
interactions.  In  preparation  for  this  discussi  n  we  r.ust  realize  that  phonons 
produced  in  the  volume  V.  are  scattered  from  this  volume  principally 
along  the  direction  q  +  qv  Therefore,  we  can  obtain  the  power  of  the 
generated  wave  simply  by  multiplying  f:q.  11  bydtw>_. 


Noncol linear  Interactions  of  Ultrasonic  Waves 

The  power  of  the  wave  Cq ^  ^  and  t*ie  power  of  the 

-> 

wave  is  (-ttui-,1.^1.  v ,N,/V) .  Using  these  expressions  and  multiplying 

fq.  11  by  fio  we  can  express  the  power  of  the  wave  generated  by  the  noncol  1  i near 
—  ►  — ► 

interaction  of  (q^,^)  and  (qn,aj^l  as 
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The  wave  of  power  P^  is  scattered  from  the  volume  V  prmcipallv  along  the 
x-direction  which  is  parallel  to  q^  +  q,„  By  examining  the  relations  l"!  and 
(N),  it  can  be  shown  that  most  of  the  phonons  (q  :  are  scattered  into  the 
solid  angle  given  by  the  first  minimum  of  the  function  2.  The  angular 


positions  of  these  minima  are  given  by  b 


(max  ) 


and  0 


(maxi 


q  .  1. 
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so  that  the  solid  angle  in  which  most  of  the  phonons  are  observed  is 
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At  large  distances  1.^,  the  average  intensity  of  the  wave  scattered 

from  the  volume  V.  can  he  calculated  as 
i 
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(14) 


where  I  and  I7  are  the  intensities  of  the  two  interacting  waves  (q^,u^)  and 
(q7,w0).  The  apparent  difference  between  (12)  and  (14)  must  be  considered 
carefully  when  experimental  data  are  taken. 

For  comparison  with  experiment  we  shall  derive  the  formula  for  the 
amplitude  of  the  ultrasonic  wave  (q  ,u_).  Using  relations  (7)  and  (Si  and 
putting  a  =  0,  a  =  q„0  1,  /2,  and  a  =  q_0  L  / 2,  the  transition  probability 

1  h  x  y  a  y  y  z  3  z  z 

rate  for  scattering  of  phonons  (q_,w_)  into  a  certain  solid  angle  AO  is  given 
by 
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where 
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liquation  (13)  gives  the  average  number  of  phonons  (q,,w,)  scattered  per  unit  ol 
time  into  “In'  solid  angle  All  along  the  direction  determined  by  the  deviation 
angles  and  Thus,  by  multiplying  llq .  (13)  by 'flu..  and  dividing  it  by 

K“V. .,  oiK’  can  ciilculate  the  average  intensity  of  the  scattered  wave  at  a 
distance  R  from  the  interaction  volume.  From  this,  one  can  obtain  the  amplitude 


of  the  ultrasonic  wave  in  the  form 
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vw = 
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(17) 


This  formula  gives  the  distribution  of  the  amplitude  of  the  generated  wave  as  a 

function  of  9  and  9  .  It  gives  an  excellent  possibility  to  analyze  experimental 
y  z 

data  in  detail  and  to  correct  for  geometrical  diffraction  of  the  generated  wave. 

It  is  interesting  to  compare  Eq.  17  with  similar  formulae  obtained  in 
references  9  and  11.  Agreement  is  obtained  if  we  put  into  Eq.  17,  o  =  1  and 
R  =  ,  where  S  is  the  cross  sectional  area  of  the  wavefront.  However,  making 

this  substitution  eliminates  the  "wave"  character  of  the  scattering,  and  thus 
leads  to  inexactness.  The  results  in  9  and  11  are  inexact  in  that  the  scattering 
process  has  been  calculated  as  though  it  were  spherically  isotropic  scattering 
with  strict  fulfillment  of  the  law  of  conservation  of  momentum.  No  such 
assumption  is  made  in  deriving  Eq.  17. 

Col  linear  Interactions:  Second  and  Third  Harmonic  Generation 
by  Sinusoidal  Ultrasonic  Waves 

Col  linear  interactions  of  ultrasonic  waves  bring  some  peculiarities  to 
the  problem  of  correct  interpretation  of  the  theoretical  formula  because 
the  created  phonons  (q.,w^)  remain  in  the  interaction  volume:  their  phase 
velocity  and  direction  of  propagation  arc  the  same  as  the  phase  velocity  and 
direction  of  propagation  of  the  interacting  waves  (we  neglect  the  effect  of 
dispersion).  First,  we  recall  that  the  optimum  conditions  for  interactions 
arc  given  by  simultaneous  fulfillment  of  the  laws  of  conservation  of  energy 
and  of  momentum.  The  corresponding  relation  can  be  written  symbolically  as 


5? 


(qj.Wj)  +  (q2,w2)  =  (qj,^) 


We  can  now  introduce  the  symbolical  vector  =  (%,a)n^’  an<*  using  it  we  can 
graphically  solve  the  relation  (18)  in  the  (q,u>)  plane  as  shown  in  Figure  1. 
Dashed  lines  in  Figure  1  should  be  drawn  so  that  the  tangent  of  the  angle  of 
deviation  from  the  q-axis  is  equal  to  the  corresponding  phase  velocity.  Thus, 
using  a  diagram  like  Figure  1,  it  can  be  shown  that  col  linear  interaction  can 
exist  when  the  interacting  waves  have  the  same  phase  velocities,  or  only  in 
special  cases  when  the  phase  velocities  are  different;  i.e.,  when  vj  ^  v2  ^  v3* 
This  relation  can  be  satisfied  only  in  directions  in  a  crystal  which  are  not 
pure-mode  directions.  Thus,  collinear  interaction  of  waves  of  different  modes 
can  exist  only  when  the  waves  are  propagating  along  a  direction  which  is  not 
a  pure-mode  direction. 

Let  us  now  consider  the  case  when  v^  =  v?  =  v^J  i.e.,  the  collinear 

interaction  of  waves  of  identical  modes  of  polarisation.  The  experimental 

situation  in  this  case  can  be  described  as  follows:  Two  ultrasonic  pulses 

(q^,Wj)  and  (q^.w^)  are  simultaneously  launched  into  the  sample  at  time  t  =  0. 

The  pulses  travel  along  a  pure  mode  direction  in  the  sample  producing  phonons 

(q^.Wj)  which  remain  in  the  space  simultaneously  occupied  by  both  pulses,  as 

we  have  assumed  v^  =  v,  =  v,.  Describing  this  situation  in  terms  of  continuous 

waves,  we  can  say  that  the  concentration  (number  per  unit  volume)  of  phonons 

(q^,w^)  at  time  t  =  L^/v.  produced  in  the  pulse  regime  will  be  the  same  as 

the  concentration  produced  by  continuous  waves  propagating  a  distance  L  . 

Using  this  conclusion  and  using  F.q.  11  for  the  concentration  of  phonons  (q_,io.) 

at  a  distance  L  ,  we  obtain 
x 
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N3(Lx>  = 
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where  v  -  v^  -  v2  -  v^.  Realizing  that^ifwNL  is  the  energy  density  of  phonons 

in  the  wave  identified  by  the  subscript  i,  one  can  calculate  the  amplitude  of 

the  wave  at  the  distance  L  as 
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a3(lx)  =—  2  qlq2AlA2Lx  (20) 
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In  deriving  both  N^L^)  and  A3(LX)  we  have  assumed  that  Nj  and  N2  are  constants 
and  that  =  0  at  t  =  0. 


Second  Harmonic  Generation 

Equations  19  and  2(1  can  be  applied  directly  to  second  harmonic  generation 
by  an  initially  sinusoidal  ultrasonic  wave.  Simply  putting  =  N2  =  NQ/2 
and  assuming  vj  =  v2  =  v3  =  v  as  before,  we  find  that  the  second  harmonic 
generated  by  a  wave  which  has  a  concentration  of  phonons  NQ  is,  from  Eq.  19: 
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where  Ns  is  the  concentration  of  phonons  in  the  second  harmonic,  while 

US  =  2u0  is  its  an8ular  frequency.  In  a  similar  manner,  from  Eq.  20  one 

finds  that  the  amplitude  of  the  second  harmonic  at  a  distance  L  from  the 

x 

beginning  of  the  sample  is 
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(22) 


This  formula  agrees  with  that  derived  from  classical  considerations.1  In  lief.  1, 

3K7+K3  2 

the  amplitude  of  the  second  harmonic  is  given  as  -  (— ^ - )A  k“a,  where  K.  and  K_ 

2  3 

t. 

arc  the  appropriate  combinations  of  second-order  and  third-order  elastic  con¬ 
stants,  respectively,  A  is  the  initial  fundamental  amplitude,  k  is  the 
propagation  constant  2tt/X,  and  a  is  the  sample  length.  Except  for  notational 
differences,  these  two  expressions  are  identical. 


Third  Harmonic  Generation 

The  third  harmonic  of  an  initially  sinusoidal  ultrasonic  wave  is  built 
up  by  two  processes:  First  is  a  three-phonon  process  in  which  the  fundamental 
phonon  (q^,^)  and  the  second  harmonic  phonon  (qs,tos)  annihilate  to  create 
the  third  harmonic  phonon  (q  ,00  ),  where  w  =  to  +  to  =  3to  ;  and  second  is  a 

t  t  u  S  O  O 

four-phonon  process  in  which  three  annihilating  fundamental  phonons  (ll0»“0) 
create  the  third  harmonic  phonon  (q  ) .  We  shall  consider  here  only  three- 
phonon  processes  which  will  give  the  significant  terms  found  in  the  third 
harmonic  amplitude.  The  four-phonon  processes  are  ignored  in  the 
following. 


Generation  of  the  third  harmonic  gradually  increases  with  the  increase 
of  the  generation  of  the  second  harmonic.  The  second  harmonic  concentration 
at  a  distance  from  the  end  of  the  sample  is  given  by  Eq.  21.  This  relation¬ 
ship  is  used  in  Eq.  ll,  along  with  the  expression  dN^/dt  =  v  dN^/dx  to  obtain 
the  concentration  of  the  third  harmonic  in  the  form 
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Applying  the  boundary  condition 
can  be  expressed  as 


1 


0  at  x  =  0,  we  find  that  the  third  harmonic 


A.  =  2A2/A 
t  s  o 


(24) 


which  is  in  very  good  agreement  with  the  formula  obtained  from  considering 
the  harmonic  generation  as  a  nonlinear  process.  To  facilitate  comparison,  we 
repeat  Eq.  1  of  Ref.  3  with  slight  notational  changes: 
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(24a) 


Comparison  of  Eqs.  24  and  24a  leads  to  the  conclusion  that  our  neglect  of  four-phonon 
processes  has  dropped  terms  corresponding  to  the  small  term  involving  fourth- 
order  elastic  constants.  This  term  was  found  to  be  negligible  in  the  experi¬ 
ments  reported  in  Ref.  3. 


IV.  PARAMETRIC  DOWN-CONVERSION  OF  TWO  ULTRASONIC  WAVES 
AND  PARAMETRIC  AMPLIFICATION  OF  AN  ULTRASONIC  WAVE 
TllROUC.il  ITS  NONLINEAR  INTERACTION  WITH  ANOTHER 
ULTRASONIC  WAVE 


The  parametric  down-conversion  mixing  process  is  observed  as  the 
generation  of  an  ultrasonic  wave  of  frequency  by  two  ultrasonic  waves  and 
oo-,  such  that  This  process  can  be  described  as  a  process  in  which 

a  phonon  (q^,w^)  is  annihilated,  and  phonons  (q7,u-,)  and  (q^w^)  are  created. 
The  process  involving  the  annihilation  of  the  phonon  (q^,w^)  and  the  creation 
of  two  new  phonons  is  suggested  in  the  Hamiltonian  density  (Eq.  3)  by  terms 
containing  3^  3*,  3*..,..  Although  spontaneous  annihilation  of  a  phonon  has  a 
very  low  probability  of  occurrence,  the  disintegration  can  be  stimulated  by  the 
presence  of  one  of  the  two  components  resulting  from  the  annihilation. 
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Let  us  assume  that  the  annihilation  of  the  phonon  is  stimulated 

by  the  presence  of  the  phonon  We  assume  that  there  are  N  phonons 

(qj,o>j)  and  phonons  (q^ ,012)  i-n  the  volume  of  normalization  V  under  con¬ 
sideration.  The  amplitude  of  the  transition  from  the  initial  state  |N  ,N-,,0> 
to  the  final  state  |N  -1,^+1, 1>  is  given  also  in  this  case  by  Eq.  4,  provided 
that  Nj  and  N.,  are  much  larger  than  unity.  Repeating  the  procedure  described 
in  the  previous  section,  one  can  again  derive  Eq.  11.  This  means  that 
parametric  down-conversion  is  described  by  the  same  formulae  as  the  parametric 
up-conversion;  however  the  down-conversion  is  accompanied  by  the  effect  of 
the  amplification  of  the  stimulating  wave.  The  amplification  of  the  stimulating 
wave  results  from  the  fact  that  the  phonons  created  from  the  disin¬ 

tegration  of  the  phonons  (q^,!^),  contribute  to  the  stimulating  wave. 

Let  us  now  calculate  the  effect  of  the  amplification,  and  simultaneously 
take  into  account  the  attenuation  of  the  stimulating  wave.  We  define  the 
attenuation  coefficient  of  an  ultrasonic  wave  as 


1  dN  1 

2N  dt  2N  dx  V  aV 


(25) 


where  a',  the  temporal  attenuation  coefficient,  expresses  the  decrease  of 
amplitude  per  unit  of  time,  and  a  expresses  the  decrease  per  unit  of  length. 
Using  this  definition  we  can  modify  Eq.  1 1  to  include  attenuation: 


rlT  C. 


V2q3  "l"2"3 
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2  2 

V1V2  3 


3  LxNl 


2«; 


(20) 


Assuming  that  N,  =  0  at  t  =  0,  the  energy  density  of  the  stimulating  wave  at 

time  t  =  L  /v,  is 
x  2 
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P2(LxJ 


p2(0)oxp 


qlq2q3w2u3  ,2 
LxP 
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5  2  3  3 
V1V2V3 


-  2a.l 


(27) 


where  is  the  length  of  the  sample  in  the  direction  of  the  wave  vector 

-V 

q3  =  ql  "  q2° 


Influence  of  Attenuation  on  Nonlinear  Interactions 
of  Ultrasonic  Waves 


We  have  tacitly  assumed  that  the  attenuation  of  the  interacting  waves  is 
in  all  previous  sections  except  for  the  section  on  parametric  down-conversion. 
This,  however,  never  is  completely  true,  and  therefore  this  section  will  be 
devoted  to  discussion  of  the  possible  effect  of  ultrasonic  attenuation  on 
nonlinear  interactions  of  ultrasonic  waves. 

As  before,  the  total  transition  probability  rate  (Eq.  li)  is  equal  to 
the  total  number  of  transitions  per  unit  of  time;  i.e.,  it  expresses  the 
average  number  of  phonons  created  in  the  volume  \\  per  unit  of  time. 

With  the  help  of  Eq.  25  and  Eq.  11  we  can  write  for  the  concentration  of 
phonons  (q,,^)* 

,2 


dN 

It 


“ft  G 


N3  qiq2q3  wlt02U3  ,  v  m  ,„-N 

=  - - f  2~2~5  x  1  2  "  "“S'  2 

8P  vlv2V3 


(28) 


If  we  assume  that  N, 


-2i  t 

N  .  e  and  N. 

U  1  *. 


~2azx 

e  ,  when  N  and  N  are 


I  ‘’01  '  “  2  02  '  01  "0 

the  concentrations  of  phonons  at  time  t  =  0,  Eq.  28  has  a  solution  of  the 


ero 


form 
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-h  tr 

q,q->q, 

N3  =  “  TT1  LxN01N02 
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Expanding  the  exponential  functions  and  keeping  only  the  linear  terms,  one  has 
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which  is  valid  for  small  attenuation  and  short  times  t.  From  this,  one  sees 

that  attenuation  can  be  neglected,  and  Eq.  30  reduces  to  Eq.  19,  when 

L 

x 

(a'  +  a'  +  a')  t  =  (a,  +  a.+  a_)  —  <<  1.  In  materials  with  a  low  value  of 

v  1  2  3  ^  1  2  3  v 

2 

attenuation,  the  attenuation  coefficient  often  is  expressed  by  a  =  A  to  ,  where 

-20  -17  -1  2 

A  is  a  constant  of  the  order  of  magnitude  10  to  10  cm  sec  .  Using 


2  2  2  v 

this  condition,  the  attenuation  is  negligible  if  Uj  +  << 


This 


condition  usually  is  very  well  fulfilled  when  noncollinear  interactions  at 
frequencies  between  10  and  100  MHz  are  investigated.  Therefore,  we  shall  go 
into  greater  detail  only  on  collinear  interactions. 

Collinear  interactions.  Since  the  power  of  the  ultrasonic  wave  can  be 

written  P.  =diw.N.L  L  v.  ,  wc  can  use  Eq.  29  to  write 
i  ]  i  y  z  i 
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Here  is  interpreted  to  mean  the  effective  interaction  length,  while  t  is  the 
effective  interaction  time.  Although  this  expression  is  limited  to  collinear 
interactions,  it  nevertheless  is  quite  general  in  that  it  is  valid  for 
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describing  the  interaction  of  ultrasonic  waves  of  different  polarization 
modes . 

If  the  interacting  waves  are  of  the  same  polarization  mode  we  can  again 
put  Vj  =  v2  =  v,  =  v  and  write 


G  n2 

q,q_q  7 

>  =  _ 1  -2.  -3— 1  L  Pn1Pn0 

■t  3  9  x  01  02 
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where  a  =  a  L  is  the  total  attenuation  of  each  wave  over  the  sample  length  I. 
n  n  x 

From  this,  the  amplitude  at  the  end  of  a  sample  of  length  is  found  to  be 
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(33) 


Stabilization  Distance 

At  a  certain  distance  of  propagation  in  the  sample,  the  harmonic  growth 

will  be  equal  to  the  attenuation,  with  the  result  that  a  stabilization  of  the 

waveform  will  occur.  This  stabilization  distance  L  can  be  calculated  from 

max 

F.q .  32  by  taking  the  derivative  and  setting  it  equal  to  zero.  For  a  sufficiently 
long  sample  the  stabilization  distance  is 

£n(a,  +  a_)  -  5-n  a, 

L  =  vt  =  ,,  1---- - .  C34) 

max  max  2(a^  +  a2  “  a3^ 

This  is  the  distance  at  which  the  ultrasonic  wave  q^  resulting  from  the 


nonlinear  interaction  has  its  maximum  value 


Comparing  this  equation  with  that  published  in  Ref.  b,  one  finds  a  differemu 
of  a  factor  of  two.  The  origin  of  this  difference  i  .  not  known  at  present. 


Second  Harmonic  Generation 

Assuming  that  =  A^7  =  Aq//2,  the  case  of  second  harmonic  generation 
by  an  originally  sinusoidal  ultrasonic  wave  of  amplitude  A^  can  be  calculated. 
From  Eq.  33  we  find 
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Iq2q3  2.2 
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(33) 


where  a  and  a  denote  total  attenuation  of  the  primary  wave  and  the  second 
os 

harmonic,  respectively. 


V.  CONCLUDING  REMARKS 

We  have  discussed  in  some  detail  the  nonlinear  effects  arising  from  the 
interaction  of  two  ultrasonic  waves.  As  has  been  mentioned,  there  is  good 
agreement  between  our  results  and  the  results  in  Refs.  1,  2,  3  and  15,  if  one 
compares  the  particular  limiting  cases  of  our  general  treatment  with  the  result 
in  these  papers.  This  agreement  is  a  valuable  result  because  it  shows  that  in 
the  correspondence  limit  the  quantum  mechanical  description  of  harmonic 
generation  gives  the  same  results  as  classical  calculations  for  both  the  second 
harmonic  amplitude  and  the  third  harmonic  amplitude.  This  emphasizes  the  point 
that  nonlinear  interactions  arc  a  unique  and  powerful  tool  for  the  measurement 
of  the  nonlinear  properties  of  materials. 


Finally,  we  shall  discuss  the  limiting  assumption  concerning  the  power 
of  the  generated  wave.  In  discussing  the  nonlinear  ultrasonic  interactions 


(v 

wo  have  sumed  th.it  tin-  iinti.il  state  of  t  ho  created  phonons  is  empty.  This 
■no. that  on  1  ■  the  boo,  1  an  1  no,  of  the  generat  i on  pr-'c  ss  has  been  considered. 
Hii'.  is  iaui  1 1  ■'  the  situation  dosori  boil  in  a  previous  paper  based  on  classical 
calculations.^  Hut  i a  the  present  situation  this  does  not  put  an  essential 
limitation  on  our  consideration-.  If  the  initial  state  of  the  created  phonons 
were  not  initially  empty,  the  "baeh  scattering"  should  be  taken  into  account. 
This  situation  can  be  described  s vmbo 1 i ca I  1 e  hv 


<N, 


(36) 


=  N]N2  -  (N1*Nr)N. 


The  situation  described  in  our  discussion  is  obtained  if  N,  =  N,  >>  N_.  It 
can  he  shown  that  the  correction  for  "back  scattering"  is  about  2°  when  the 
power  of  the  generated  wave  is  about  1%  of  the  primary  waves.  In  most 
practical  cases  the  power  of  the  generated  wave  is  much  lower  than  1"  of  the 
power  of  the  primary  wave.  Therefore,  the  equations  we  have  derived  arc 
adequate  to  describe  most  experimental  situations. 
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relation  (2)  are  stiffened,  and  also  there  is  a  contribution  to  the  stiffening 
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